Theory of double resonance magnetometers based on atomic alignment 



Antoine Wei^ and Georg Bison 
Physics Department, University of Fribourg, Chemin du Musee 3, 1700 Fribourg, Switzerland 



O 
O 

(N 

Oh! 



43 



o 

c3 
o 

>Y 
43 

(N 
> 

cn : 

in ■ 
o ■ 

: 

o . 

o : 

43 : 

Ph. 



X 



Anatoly S. Pazgalev 

Physics Department, University of Fribourg, Chemin du Musee 3, 1700 Fribourg, Switzerland and 
Ioffe Physical Technical Institute, Russ. Acad. Sc., St. Petersburg, 194021, Russia 

(Dated: May 26, 2006) 

We present a theoretical study of the spectra produced by optical/radio- frequency double reso- 
nance devices, in which resonant linearly polarized light is used in the optical pumping and detection 
processes. We extend previous work by presenting algebraic results which are valid for atomic states 
with arbitrary angular momenta, arbitrary rf intensities, and arbitrary geometries. The only restric- 
tion made is the assumption of low light intensity. The results are discussed in view of their use in 
optical magnetometers. 
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I. INTRODUCTION 

Since the 1950s the combination of resonant optical ex- 
citation and magnetic resonance has been an extremely 
valuable tool for atomic spectroscopy. This double res- 
onance technique Q] has not only proven useful for in- 
vestigating atomic structure, for measuring properties of 
atoms, their constituents, and their interactions, but has 
also led to important applications in atom cooling, op- 
tically pumped frequency standards, and optical magne- 
tometers. 

Magnetometers based on double resonance in atomic 
samples measure the modulus of an externally applied 
magnetic field Bo via the Larmor precession frequency 
of the sample's magnetization in that field 0, El • The 
sample is typically a vapor of paramagnetic atoms (or 
diamagnetic atoms excited to a metastable state with 
an orbital angular momentum) sealed in a glass cell. A 
macroscopic magnetization is created in the vapor by op- 
tical pumping with polarized resonance radiation. The 
magnetization precesses in the magnetic field Bq to be 
measured (referred to as the offset field) and that preces- 
sion is driven by a (much weaker) magnetic field Bi(t) 
(referred to as the rf field), co-rotating with the magne- 
tization around the offset field. 

Since the optical properties of the medium, character- 
ized by its complex index of refraction, depend on its 
spin polarization, the driven magnetization will induce 
periodic modulations of those properties [J, |5j , which are 
then detected. In most applications, the same light beam 
used to polarize the medium is also used to monitor the 
oscillations by measuring either the power or the polar- 
ization state of the transmitted beam. The frequency 
of the induced oscillations coincides with the oscillation 
frequency u> of the rf field, and their amplitude depends 
in a resonant way on the detuning between u> and the 



Larmor frequency luq — ^fBq associated with the offset 
field. The Lande factor jp — gpfis/fiis characteristic for 
the pumped atomic state with total angular momentum 
F. 

Most practical double resonance devices rely on atomic 
orientation prepared by optical pumping with circularly 
polarized light. In this paper we present a theoretical 
study of the resonance signals obtained in double res- 
onance spectroscopy using linearly polarized light. As 
shown first by Bell and Bloom [(J magnetic resonance in 
aligned media leads to signal modulations at the fun- 
damental and at the second harmonic of the rf fre- 
quency. We derive algebraic expressions for the spectral 
line shapes of the in-phase and quadrature components 
of both signals and their orientation dependence. Pre- 
vious theoretical treatments of such signals [21 |8j were 
restricted to specific angular momentum states (J = 1) 
or to low rf powers. The results presented here are more 
general in the sense that they apply to arbitrary spin sys- 
tems and that they are valid for arbitrary rf power levels 
and for arbitrary orientations of Bq with respect to the 
light polarization. 



II. POLARIZED ATOMIC MEDIA 

A. Atomic Multipole Moments 

The density matrix p of an ensemble of polarized atoms 
with angular momentum F can be expressed in terms of 
atomic multipole moments m^ q according to |j| 
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where the Tq are standard irreducible tensor operators 

F F 

^ k) (F) = £ £ {-i) F - M ' (2) 

M=-F M' = -F 

x(F,M,F,-M'\k,q) \F,M) (F,M'\ 

constructed from the angular momentum states \F,M), 
and where the multipole moments rrikq arc defined by 

m k>g = (T^)=Tr(pT^). (3) 

The three multipole moments TOfc=i,g=-i,o.+i represent 
the orientation of the medium, while the five components 
TOfc=2,g=-2,...,+2 represent its alignment. The multipole 
moments nik,q=o are called longitudinal multipole mo- 
ments and their value depends only on sublevel popula- 
tions. The multipole moments mk,q^o represent sublevel 
coherences and are called transverse moments. The rep- 
resentation of the atomic polarization in terms of mul- 
tipole moments has a significant advantage over a rep- 
resentation in terms of sublevel populations and coher- 
ences. In principle, both representations require the same 
number of parameters for the complete description of the 
atomic ensemble. However, because electric dipole radi- 
ation couples only jToj to orientation (fc = 1) and align- 
ment (k = 2) it is sufficient to specify the corresponding 
3 + 5 = 8 multipole moments for the complete descrip- 
tion of the system's optical properties. Moreover, specific 
light polarizations couple only to specific subsets of these 
8 multipole moments, so that the use of the tensor for- 
malism in systems with large angular momenta leads to a 
significant simplification of the mathematical treatment. 
In the case discussed here only the (real) multipole mo- 
ment m2,o wm be relevant. This approach therefore al- 
lows one to derive results valid for systems with arbitrary 
angular momenta. 

A resonant circularly polarized laser beam interacting 
with an unpolarizcd atomic sample will create orientation 
(k = I, vector polarization) and alignment (fc = 2, tensor 
polarization) in the sample by optical pumping. The low- 
est order multipole that a linearly polarized light field can 
create is an atomic alignment. While only atomic states 
with J > 1/2 can be oriented, the condition J > 1 has 
to be fulfilled for the creation of an aligned state. Note 
that an alignment along the direction of light propaga- 
tion can also be produced by pumping with unpolarized 
light 0, The ground state of alkali atoms has an 

electronic angular momentum J = 1/2, which cannot 
be aligned. However, the hyperfine interaction with the 
nuclear spin splits the ground state into two hyperfine 
levels with total angular momenta F± — I ± J, which 
can be aligned provided F > 1. An alignment can there- 
fore be prepared and/or detected only if the light source 
has a sufficient spectral resolution to excite a single hy- 
perfine transition. In general the Dopplcr (and pressure) 
broadened spectra of discharge lamps, used in conven- 
tional optically pumped magnetometers (OPM), cannot 
be used to address individual hyperfine lines and hence 




FIG. 1: Parametrization of the DRAM geometry, in which 
Bi is perpendicular to Bo. The rotating wave (rw) frame and 
the rotating magnetic field Bi are shown for one particular 
moment in time (t = 0). 



do not allow one to create nor to detect a ground state 
alignment. However, radiation from a narrowband laser 
can resolve the hyperfine structure and it is well known 
that a linearly polarized laser beam can create an atomic 
alignment. 



B. DROMs and DRAMs 

We will refer to an OPM based on atomic orientation 
as DROM (double resonance orientation magnetometer) , 
while we will speak of DRAM (double resonance align- 
ment magnetometer) when the magnetization has the 
symmetry of an atomic alignment. Most of the past re- 
search work on double resonance spectroscopy dealt with 
oriented vapors, although alignment induced by (unpo- 
larized) lamp pumping in J = 1 metastable states of He 
was already reported in 1960/61 |6J,lllJ. In alkali atoms, 
alignment produced by lamp pumping can be observed 
using line splitting by the quadratic Zeeman effect |T^ or 
isotope filtering yjIT^]- The latter technique is, however, 
restricted to Rb and cannot be applied to other alkalis. 
As mentioned above, linearly polarized laser radiation is 
an efficient means for producing alignment and a discus- 
sion of linearly polarized laser pumping in metastable 4 He 
can be found in 0, El ■ These authors have investigated 
several magnetometry techniques using both orientation 
and alignment signals and they observed magnetic res- 
onances involving alignment using rf fields, light inten- 
sity modulation, polarization modulation, and frequency 
modulation. A variant of the latter technique - in which 
the transmitted lights' polarization, instead of intensity, 
was measured - was realized with 87 Rb 0, 0] . 



III. MODEL CALCULATIONS 

A. DRAM geometry 

We will restrict the discussion to geometries in which 
the rf field B\ is perpendicular to the offset field Bq. Be- 
cause the signals are independent of the orientation k, the 
geometry of the problem is fully determined by Bq, Bi, 
and <Tso that one can consider the parametrization shown 
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in Fig. ^ in which 7 denotes the angle between the off- 
set field and the light polarization [8|]. The orientation 
dependence of the signal amplitudes is given by 7, while 
the polar angle a will lead to a phase shift in the time 
dependence of the oscillating signals (Sec. HVB(I . 

In Fig. ^ we have applied the rotating wave approx- 
imation by decomposing the rf field into two counter- 
rotating fields, from which we have retained only the 
component co-rotating with the alignment. This com- 
ponent is shown in its position at time t — 0, thereby 
defining the time origin of the phases of the oscillating 
signal components. The magnetometer signals are cal- 
culated following the 3-step (prepare, evolve, probe) ap- 
proach introduced in 0, 0j and discussed in detail in 
. In the first step of this model (preparation) one as- 
sumes the existence of a given alignment in the atomic 
medium, without specifying how this alignment was cre- 
ated. Details of the preparation process (optical pump- 
ing, collisions with electrons or ions, spin exchange, etc) 
thus do not need to be known. In the second step (mag- 
netic resonance) the initial alignment is allowed to evolve 
towards a steady state value determined by the interac- 
tion with the external fields and relaxation. Finally one 
considers in the third step (probing) how the steady state 
alignment affects a linearly polarized light beam travers- 
ing the medium. Strictly speaking this approach is only 
valid for pump-probe experiments, in which the atoms in- 
teract with spatially or temporally separated light fields 
and where the equilibrium of step 2 is reached "in the 
dark". However, as shown previously |l6l Il7j| for level 
crossing signals the results obtained from the 3-step ap- 
proach give an excellent description of experimental find- 
ings if the light intensity is sufficiently weak. Limitations 
of the model will be addressed in section IIII El 



B. Step 1: Alignment creation 

We describe the alignment created by the preparation 
process in a coordinate frame where the quantization axis 
lies along the light polarization (lab frame in Fig.^l. In 
that frame the only non-vanishing alignment component 
created by optical pumping with linearly polarized light 
is the longitudinal multipole moment mifp which can be 
expressed in terms of the sublevel populations pu as 



mZ = N 2 (F) 
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where N%(F) is a normalization constant [9j. In the pres- 
ence of an offset field Bq the alignment components per- 
pendicular to Bq will relax to zero yielding a steady state 
value of 
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for the alignment along the magnetic field, given by the 
projection of m^p on the field direction. Note that this 



steady state is reached only when the Larmor frequency 
is much larger than the transverse relaxation rates. This 
condition is well fulfilled for high-Q magnetic resonances 
while it is not met by zero-field level-crossing resonances 
(ground state Hanle effect, nonlinear Faraday effect [lql- 



C. Step 2: Magnetic resonance 

This step describes the magnetic resonance process, 
i.e., the evolution of the alignment under the combined 
actions of the magnetic fields Bq, B±, and relaxation 
processes. It is described in a coordinate frame, which 
is related to the lab frame by a static rotation of —7 
around the y-axis and then by a dynamic rotation, at 
the frequency u, around the new z-axis. In this frame, 
which we call the rotating wave frame (rw frame), the 
offset field Bo is along z, while the rf field £?i appears 
to be static and oriented at an angle a with respect to 
the x-direction (see Fig^l. This is the usual field con- 
figuration for describing magnetic resonance processes. 
Note that m^Q is not affected by the transformation 
to the rw frame. Due to the rotation of the coordi- 
nate frame, a fictitious magnetic field B f = —lJz/jf 
appears in the rw frame, and the atoms see a total field 
-Btot = B\ cosax + B\ sinay+ (Bq — uj/j)z. 

The evolution of the the system's density matrix is 
described by the Liouville equation 



dt 



(5) 



with H = —fi - Btot 1 and where p r ciax describes the relax- 
ation processes. Inserting the multipole decomposition 
||TJ into JSJ) yields the following equations of motion for 
the multipole moments rni, q 



dt 
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in which u>x — jfBi is the Rabi frequency of the rf field, 
p± = exp(±ict) are phase factors that describe the ori- 
entation of u>i in the xy-plane, and S = u> — u>q is the 
detuning of the radio frequency u> with respect to the 
Larmor frequency loq. The relaxation terms are given by 
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where we assume that the multipole moments m 2 , g re- 
lax at rates r| g |. The transverse alignment components 
m 2jq ^o relax towards a zero value, while the longitudinal 
component to 2 q relaxes towards rn 2 Q , introduced in step 
1. Below, we present very general results for the case 
when all three relaxation rates are different, although it 
is well known that there are specific relations between 
those rates when the relaxation mechanism has specific 
rotational symmetries 0] . Eqations are a generaliza- 
tion of the well known Bloch equations 



q' 



relax 
L l,q 



-1,0,1, (8) 



describing the evolution of the three orientation compo- 



nents TTi\q. 

We have used a computer algebra software to de- 
termine algebraic expressions for the steady state 
(d/dt m 2 , q = 0) solutions m 2>g of ©■ These solutions 
are then transformed back to the lab frame, by first ap- 
plying a dynamic rotation at the rate — w around the 
z-axis, and then a static rotation by 7 around the y-axis 
of the rw frame. In this way one can derive algebraic 
expressions for the time dependent multipole moments 
m2,q{t) in the lab frame. 



D. Step 3: Alignment detection 



proportional to m 2 fi(t). We define the time dependent 
DRAM signal, normalized to the longitudinal alignment 
initially produced by the optical pumping, as 



S(t) 



m 2 ,a(t) 



, 1111 
l 2,0 



(10) 



E. Validity of the three step approach 

The three step approach is only valid if steady state 
conditions are reached in steps 1 and 2. This is fulfilled 
when the pump rate T p , at which alignment components 
are modified by the interaction with the light is negligible 
compared to the relaxation rates Ti q i . This condition can 
be realized experimentally at low light powers, however, 
at the cost of a decreased signal to noise ratio. OPMs 
are known to perform best when T p is comparable to IVi . 
The lowest order correction, taking the depolarization of 
light interactions into account, can be described by 



r 



r 



r 



vPl, 



11) 



where Pl is the laser power. 

It is well known that substitution (|llfl is valid to all 
orders in Pl for a DROM in a spin 1/2 system [2J}, in 
which orientation is the only multipole moment that can 
be created. It is reasonable to assume that the same 
statement can be made for a DRAM in a spin 1 system 
in which alignment is the only multipole moment created 
and detected by the linearly polarized light. For angular 
momenta F > 1 the creation of higher order (k > 2) mul- 
tipole moments and their transfer back to (detectable) 
k = 2 moments limits the validity of substitution 1111) to 
low light powers. 



In the third and final step, we calculate the effect the 
time dependent multipole moments have on the optical 
absorption coefficient of the medium. One can show that 
the absorption coefficient of a medium described by mk, q 
for linearly polarized light is proportional to 

A) /2 , 

K (X ~y| TO 0,0 - y 3^2^2,0 , (9) 

where the (analyzing power) Ak depends only on the 
states \n g , L g , J g , F g ) and |n e , L e , J e , F e ) coupled by the 
light. The multipole moments in are defined with re- 
spect to a quantization axis oriented along the incident 
light polarization, which is the case in the lab frame, i.e., 
the frame in which the results of step 2 are expressed. 
The monopole moment mo,o describes the total popula- 
tion of the hyperfine ground state \n g ,L g ,J g ,F g }. We 
assume the optical transition to be closed and the light 
intensity to be so weak that excited state populations 
remain negligible, so that the monopole moment does 
not depend on time. The only time dependent (oscillat- 
ing) component of the absorption coefficient is therefore 



IV. RESULTS 

The calculation outlined above yields signals S^, (t) and 
S2U1 (t) which are modulated at the rf frequency u> and at 
its second harmonic 2w, and which can be written as 

S u (t)= M7)[ D u cos (wi - a) (12a) 

— A w sin (cut — a)] , 
S 2 u (t) = h 2ul (7) [ -Mu cos (2wt - 2a) (12b) 
-Diu, sin (2ut - 2a)] , 

where the angular dependence of the signals ^(7) and 
^(7) (Fig-0) is given by 

3 

^(7) = - sin 7 cos 7 (3 cos 2 7 — l) (13a) 
hauii) = ^ sin 2 7 (1-3 cos 2 7) . (13b) 

As stated earlier the orientation angle a of the rf field 
appears as a phase shift. The in- phase and quadrature 
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FIG. 2: Angular dependence of the first (hu, left) and second (hzu, right) harmonic signals on the angle 7 between the light 
polarization e and the magnetic offset field Bo- 



components of the signal have both absorptive {A Ul A 2u ) 
and dispersive (D^, D 2uj ) lineshapes given by 
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r (r 2 + 5 2 )(r 2 
+ (r^ (2r + 3r 2 ) - as 2 (r - 3r x )) w 2 
+ (r + 3r 2 K, (15) 



Equations (|14(l and l|15() can be simplified substan- 
tially if we assume an isotropic relaxation by setting 
To = Ti = T2 = r. We will stick to this assumption 
in the following discussion since it does not change the 
general properties of the spectra. We further simplify 
the notation by introducing a dimensionless rf satura- 
tion parameter S T i = (u>i/T) 2 and a normalized detuning 
x = 5/T. With these assumptions and definitions we 
obtain 



D u (x,S I f) 
A LU (x, S T f) 
D 2w {x,S ri ) 
A 2oJ (x,S T {) 



x(l- 2S lf + Ax 2 ) ^/B7 t 
(1 + S I{ + x 2 ) (1 + AS I{ + Ax 2 ) 
(1 + S lf + Ax 2 ) ysvT 

(1 + Srf +X 2 )(l+ASri +Ax 2 ) 
(1 + S ri + X 2 ) (1 + ASrf + AX 2 ) 

(1 + grf - 2x 2 ) S rf 
(l + 5 rf + x 2 )(l + AS l{ + Ax 2 ) 



(16a) 
(16b) 
(16c) 

(16d) 



In the following we will discuss in detail the different 
properties of the signals l|12|) with lineshapes (|16|l . 



A. Line shapes 

For low rf intensities, i.e., for SVf — > 0, the lineshapes 
(|TFj|| reduce to 



and 



D u (x,Srt -f 0) = 

A u (x,S r{ ->0) = 
DaufaSrt -> 0) = 

-Afcjfo Srf -> 0) = 



a; 

1 + a; 2 
1 

1 + x 2 



Srf , 

/ Srf > 



3a; 



(1 + a; 2 ) (1 + 4x 2 ) 

l~2z 2 
(1 +x 2 ) (1 + 4x 2 ) 



(17a) 
(17b) 

SW, (18a) 
S rf . (18b) 



Expressions l|13ll and (|17|l , correspond to results obtained 
in earlier work [7J, |8j . 

The corresponding spectra can be seen in the leftmost 
columns of Figs. |31 and 0] These figures also show how 
the lineshapes of the absorptive and dispersive signals 
Auj, D^, A 2uJ , and D 2uJ change with increasing rf intensity. 
A narrow additional spectral feature appears in the first 
harmonic signal for S r f > 0.5. The origin of this structure 
can be explained as follows: the basic interaction of the rf 
field with the angular momentum is the coupling of ad- 
jacent Zeeman sublevels. The corresponding AM = ±1 
coherences oscillate at the rf frequency u) and their de- 
tection by the light field constitutes the first harmonic 
signal. With increasing rf power, a further interaction of 
a AM = 1 coherence with the rf field B\ becomes pos- 
sible and leads to the creation of a AM = 2 coherence, 
whose oscillation produces the second harmonic signal. 
An additional interaction of the AM = 2 coherence with 
the rf field produces both AM = 3 and AM = 1 co- 
herences. While the former cannot be detected optically 
since linearly polarized light couples at most to AM = 2 
coherences, the latter directly contributs to the first har- 
monic signal. In this sense, the additional features can be 
understood as resulting from the creation of a AM = 2 
coherence by a second order interaction with the rf field, 



6 




FIG. 3: (Color online) Top line: Line shapes of the absorptive (A u ) and dispersive (D^) components of the first harmonic 
signal for different values of the rf saturation parameter S T {. Bottom line: Shape of the phase signal ip u for the same values of 
Srf. x = (u — Uo)/T is the normalized detuning of the rf frequency ui with respect to the Larmor frequency uio- 



S rf =0.3 S r ,=1 S rf =3 S rt =10 




-4 -2 2 4 -4 -2 2 4 -4 -2 2 4 -4 -2 2 4 

detuning x 



FIG. 4: (Color online) Top line: Line shapes of the absorptive (A2u) and dispersive (Z?2w) components of the second harmonic 
signal for different values of the rf saturation parameter SW. Bottom line: Shape of the phase signal ip2u> for the same values 
ofSrf. 



the evolution of that coherence in the offset field, and 
back-transfer to a AM = 1 coherence by an additional 
interaction with the rf field followed by detection of that 
coherence at the first harmonic frequency uj. 

The additional narrow feature in the absorptive signal 
that appears at large rf intensities was already observed 
in the pioneering work by Colegrove and Franken on op- 
tically induced alignment Based on the results pre- 
sented above the feature can be explained in terms of line 
superpositions. In fact, the expression for the absorptive 
first harmonic signal A^ (Eq. I16b|) can be rewritten as 



A^ — 



1 + grf 
1 + 5 rf + x 2 



Srf- 



4S rf 



1 + 4SW + Ax 2 



(19) 



i.e., as a superposition of two absorptive Lorentzian line 
shapes whose widths, for 5 r f — > 0, differ by a factor of 



2, while they become equal for large values of S t f- The 
appearance of the central dip is a consequence of the 
different rates at which the contributions saturate with 
increasing SVf. At low rf powers, the amplitudes of the 
two contributions to the first harmonic signals (|19|l grow 

as Sl( 2 oc u>% and SJ 3 / 2 oc u>f respectively, which reflects 
that these resonances correspond to first and third or- 
der processes as discussed above. The second harmonic 
signals (Eas. ITSal I18b(l . on the other hand, grow as u 2 , 
which reflects their second order nature. 

The dependence of the line widths, i.e., the frequency 
separation Awpw of the maxima and minima of the dis- 
persive signals D u and Z?2w on SVf can be inferred from 
the derivatives of l|It)a|) and Ultjcfl . The dispersive line 
shape of the second harmonic signal is thus 2.6 times 
narrower than the corresponding first harmonic signal. 
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B. The phase of the signals 

We define the phases ipu, {<p2u) of the first and second 
harmonic signals as the phase difference between S u (t) 
(S2u(t)) and oscillations that are proportional to coswi 
(cos 2cut) respectively. This definition is motivated by the 
fact that the rf field in the lab frame is proportional to 
cos cot. As an alternative to the parametrization of the 
signals in terms of in-phase and quadrature components 
fEa. ll2fl one can write S u (t) and S2u (t) m terms of mod- 
uli R u , R,2 U and phases ip u , ip2 U according to 

Su, = h U) ('j)R u (x,Stf)coa[urt + ipu)(x,S T {)] , (20a) 
S 2 u = h 2 u('Y)R2u(x,S l{ )cos[2u;t + tp2u,(x,S li )] , (20b) 

with 



<f2u 



arctan 



= arctan 



-^-Lu{x, Srf ) 

D U {X, Srf) 

D 2 uj(x, 5rf) 
A 2u (x,S I f) 



2a. 



(21) 
(22) 



A dual-phase lock-in amplifier can be used to extract 
the tp's from the signal. Such amplifiers usually allow one 
to apply an offset phase ip os to the signal which then adds 
to a. For stabilization purposes, it is practical to choose 
that offset phase such that the total phase vanishes at 
the center of the resonance ip u (x = 0) = Lp2u>{x = 0) = 
0. This choice avoids phase discontinuities and provides 
a signal <f(x) oc x (for 1 < 1) that is proportional to 
magnetic field changes near the center of the resonance. 
In the parametrization defined above this can be realized 
simultaneously for both signals when a + (p os — n/2. In 
that case the phase signals take the form 



tfiuj — — arctan [ x 
ip2ui = — arctan 



1 - 2S rf + Ax 2 
1 + S rf + Ax 2 

1 + Srf - 2x 2 



(23a) 
(23b) 



Examples of tp w (x) and (x) for different values of S r f 
are shown in Figs. and 0] respectively. 

The phase of the first harmonic signal <p u (Eq. I23afl fol- 
low rf power (S r f <C 1/2) is identical to the phase of an 
orientation based magnetic resonance signal. Conversely 
to that DROM phase, which does not depend on S r f , the 
DRAM phases up^ and ip2u depend on S r f, and ip^ even 
changes the sign of its slope at Srf = 1/2. Note that (fi2w 
makes a total phase swing of 2-7T across the resonance, 
while ip^ swings only by 7r. 



V. MAGNETOMETRY CONSIDERATIONS 

There are two modes of OPM operation, the self- 
oscillating mode and the free running mode, which differ 
in the way the magnetic field value is extracted from the 
signals described above. Self-oscillating (or phase stabi- 
lized) magnetometers use feedback to keep the driving 



frequency u> equal to u>o, and in such magnetometers, 
uj is measured with a frequency counter [211]. We will 
not discuss the self-oscillating mode in detail since the 
feedback is difficult to describe analytically. Instead, we 
discuss the free-running mode which yields an identical 
magnetic field sensitivity. The goal of the following dis- 
cussion is not the derivation of absolute field sensitivities, 
but rather a comparison of the relative magnetometric 
sensitivities that one can expect from the first and sec- 
ond harmonic signals. 

As a free-running magnetometer application we con- 
sider the recording of signal changes induced by vari- 
ations of the offset field Bq while the rf frequency to 
is kept constant. Either the dispersive signals D u and 
£>2w (Eqs. I16alll6cfl or the phase signals <p u and (f2u 
(Eqs. I23all23b|) can be used as discriminating signals as 
they both feature a linear dependence on Bq changes near 
the center of the resonance (x = 0). 

The resolution with which field changes can be de- 
tected is limited by noise processes such as photon shot- 
noise, electron shotnoise (in a photodiode), and spin pro- 
jection noise, all of which have a white noise spectrum. 
We specify the magnetometric sensitivity in terms of the 
noise equivalent magnetic field (NEM) [2^, [2^], which 
is the amplitude Ai? NEM of field fluctuations which in- 
duce fluctuations of the (dispersive) signal S(Bo) that 
are equal to the signal noise AS: 



AB 



NEM 



AS (dS/dB \ x=o y 



One can show that for a given noise level of the modu- 
lated signal one obtains the same NEM either from the 
demodulated signal S = D u (D2u) or from the demodu- 
lated signal S = </? w (v2w)- The minimal value of A_B NEM 
is thus obtained under conditions which maximize h u ■ s u 
and h2uj ■ S2u , where the (on resonance) slopes s w and S2u 
are given by 



dD u 



$2* 



dx 
dD 2u} 



dx 



(1 — 2S r f)v / Srf 
x=0 ~ (l + Srf)(l + 4S rf ) 
3S r f 

x=0 ~ (l + S rf )(l + 4S rf ) 



(24a) 
(24b) 



Figure [S] shows their dependence on the rf intensity. The 
zero crossing of the slope of the first harmonic signal at 
S r f = 0.5 marks the emergence of the narrow central 
feature in Fig. 0] 

The maximal sensitivity (minimal AB^^f ) is achieved 
by choosing a geometry which which maximizes flu and 
h<2u> and an rf intensity, which maximizes s w and S2u> ■ For 
the first harmonic signal one finds max[/i w • s u ] = 0.141 
for 7 = 25.5 degrees and S r f = 0.079, while for the second 
harmonic signal one has max[ft,2aj ■ S2w] = 0.25 for 7 = 
90 degrees and S r f = 0.5. Under optimized geometrical 
and rf power conditions and for a given noise level the 
second harmonic signal is thus expected to yield a 1.8 
times higher sensitivity to magnetic field changes than 
the first harmonic signal. 
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VI. SUMMARY AND CONCLUSION 

We have presented a general theoretical framework for 
the calculation of optical rf double resonance signals that 
can be easily applied to oriented or aligned atomic media. 
The theory yields analytical results with a broad range 
of validity, and is only limited by the assumption of low 
light power. 

DRAMs, i.e., magnetometers that use linearly polar- 
ized light, present several potential advantages over the 
well known DROM scheme: 

• Line widths: The line shapes of the second har- 
monic DRAM signal have significantly narrower 
linewidths than the DROM signal under identical 
conditions. Narrow linewidths potentially increase 
the magnetometric sensitivity and suppress system- 
atic errors in optical magnetometers due to long 
term baseline drifts p4|. 

• Light shift: In DROM devices the interaction of the 
atoms with the circularly polarized laser light leads 
to M dependent energy shifts (vector light shift) of 
the sublevels \F, M) when the laser frequency is 
not centered on the optical resonance line. In that 
case, power and frequency fluctuations of the laser 
mimic magnetic field fluctuations, thereby limiting 
the magnetometric performance. In the DRAM de- 
vice the linearly polarized light produces a tensor 
light shift |25[ depending on M 2 , which does not 
have the characteristics of the Zeeman interaction 
and will therefore not affect the magnetometric per- 
formance. 

• Geometry: The DROM scheme achieves a maxi- 
mal sensitivity for 9 = tt/4. The 45 degree angle 
that the laser beam has to make with the mag- 
netic field seriously limits applications which call 
for a compact arrangement of multiple sensors. In 
multi-channel devices, as required, e.g., for cardio- 
magnetic measurements [2(J , the use of the DRAM 
signals offers the advantage that the laser beam can 



be oriented either parallel or perpendicular to the 
offset field. 

• Vector magnetometry: Both the DROM and the 
DRAM devices are scalar magnetometers and their 
resonance frequency measures |-Bo|- However, the 
DRAM signals can be used to realize a vector mag- 
netometer, since the ratio of Rzu/Rui is propor- 
tional to tan 7, for all values of S T f. The knowl- 
edge of | Bo | and 7 locates Bq on a cone and the 
variation of the signal with 4>, obtained by rotating 
the polarizer, will determine the polar angle of Bq 
on the cone. In this way the DRAM scheme can 
be used to infer all three vector components of the 
field. 

• Relaxation: We are in the process of performing 
extensive experimental studies of the DRAM prop- 
erties in paraffine coated cesium cells j^. First 
results indicate that a description of the signals us- 
ing three independent relaxation rates fEas. IT4*|) is 
required to describe the experimental lineshapes in 
detail. The DRAM signals thus seem to offer a con- 
venient way for studying spin relaxation processes 
in aligned media. 

• Signal noise: Diode lasers are convenient light 
sources for double resonance experiments. How- 
ever, they often feature a 1// (flicker) intensity 
noise at low frequencies that turns into the white 
shot noise level at higher frequencies. The detec- 
tion of Z?2uj and y4.2 W at twice the Larmor frequency 
makes it easier to operate in a region where the 
laser noise is less affected by flicker noise. 

Cs OPMs in the DROM geometry have a shot noise 
limited sensitivity of 10 fT/Hz 1 / 2 |2g. A direct exper- 
imental comparison of DRAM and DROM magnetome- 
ters is currently underway in our laboratory. This study 
will show if the potential adventages can be realized in 
practice. 
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